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Abstract 

In the last few years, new observations by CHANDRA and XMM have shown that 
Pulsar Wind Nebulae present a complex but similar inner feature, with the presence 
of axisymmetric rings and jets, which is generally referred as jet-torus structure. Due 
to the rapid growth in accuracy and robustness of numerical schemes for relativistic 
fluid-dynamics, it is now possible to model the flow and magnetic structure of the 
relativistic plasma responsible for the emission. Recent results have clarified how 
the jet and rings are formed, suggesting that the morphology is strongly related to 
the wind properties, so that, in principle, it is possible to infer the conditions in the 
unshocked wind from the nebular emission. I will review here the current status in 
the modeling of Pulsar Wind Nebulae, and, in particular, how numerical simulations 
have increased our understanding of the flow structure, observed emission, polar- 
ization and spectral properties. I will also point to possible future developments of 
the present models. 

Key words: star: pulsar, magnetohydro dynamics, relativity, ISM: supernova 
remnants, 



1 Introduction 



Pulsar Wind Nebulae (PWNe) are bubbles of relativistic particles and mag- 
netic field produced by the interaction of the ultra-relativistic pulsar wind 
with the ambient medium (ISM or SNR). The best example of PWN and the 
prototype of the entir e class is the Crab N ebula. The basic model of a PWN 
was first presented by iRees &: Gunnl (11974 ). and deve loped in more detail by 
Kennel &: Coroniti ( 1984a , bl) (KC84 hereafter), and Emmering fc Chevalier 



fll987h . The ultra-relativistic wind produced by the pulsar is confined inside 
the SNR, and slowed down to non relativistic speeds in a strong termination 
shock (TS). At the shock the toroidal magnetic field of the wind is compressed, 
the plasma is heated and particles are accelerated to high energy. A bubble 
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of high energy particles and magnetic field is produced where the post-shock 
flow expands at a non relativistic speed toward the edge of the nebula. This 
model was particularly successful and was able to explain many of the observed 
properties of the Crab Nebula. The continuous, non-thermal, very broad-band 
spectrum, extending from Radio to X-ray, with spectral index in the range 
0-1.2, steepening with increasing frequencies, was modeled as synchrotron 
emission, associated with the g yration of accelerated h i gh energy particles in 



the c o mpressed magne t ic field (IVeron-Cettv fe Woltierl. 11 993; Ba ndiera et al. 



19991 : IWeisskopf et all . l2000l ; IWillingale et all l200ll ; Mori et all 120041 ) . The 



presence of an under-luminous region, centered on the location of the pul- 
sar, was interpreted a s eviden c e for the ultra- r elatiy i stic unshocked wind - 



Polarization measures (Wilson, 1972; Velusamy, 1985|; ISchmidt et all 11979 



Hickson fe van den Berghl . Il990l ; iMichel et all Il99ll ) show that emission is 
highly polarized and the nebular magnetic field is mostly toroidal, as one 
would expect from the compression of the pulsar wind. Interestingly this 
model clearly explains why Crab Nebula appears bigger at smaller frequencie s 
flVeron-Cettv fe Woltierl . Il993l : iBietenholz et all 1 19971 : iBandiera etldl Il998h . 
High energy X-ray emitting particles have a short lifetime for synchrotron 
losses, and they are not advected far away from the TS; in contrast, the syn- 
chrotron lifetime for radio particles is longer than the age of the nebula, so 
they fill the entire volume. Based on this simple model it is possible to con- 
strain some of the properties of the pulsar wind, at least at the distance of the 
TS. To explain the dynamics of the plasma, as well as the emission properties 
of the nebula, the Lorentz factor of the wind is estimated to be ~ 10 6 , and 
the ratio between Poynting flux and kinetic energy o ~ 0.003. The model 
also explains why many PWNe (i.e. Crab Nebula, 3C58) have an elongated 
axisymmetric shape: the toroidal magnetic field in the nebula exerts a higher 
pressure on the confi ning ambient medium along the polar directi on than on 
the equatorial plane (IBegelman fe Lil . Il992l ; Ivan der Swaluwl . 120031 ) . 



In the model by KC84 the SNR has only a passive role providing the confine- 
ment of the PWN. However, the details of the flow structure in the PWN and, 
as a consequence, of the emission properties depend critically on the boundary 
condition with the SNR. Given the complexity of the interaction and the ex- 
istences of different timescales, a detailed study of the evolution of the system 
PWN-SNR, has been possib le only recently, thanks mostly to the improvement 



in computational resources (Ivan der Swaluw et all 120011 : iBlondin et all 12001 



Bucciantini et all 120031 : Ivan der Swaluw et all |2004j) . Comparing the energy 



in the SNR (~ 10 ergs) to the total energy injected by the pulsar in the PWN 
(~ 10 49 ergs), it is immediately evident that the PWN cannot significantly af- 
fect the global evolution of the SNR, but on the other hand, the evolution of 
the SNR can have important consequences for the PWN. Three main phases 
can be identified in the evolution of a PW N inside a SNR (for a more com- 
plete discussion of PWN -SNR evolution see lReynolds fe Chevalierl (119841 ) and 
Gaensler fe Sland (120061 )). In the first phase the PWN expands inside the cold 
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Fig. 1. X-ray CHANDRA images (NASA/CXC/SAO) of the inner region of PWNe: 
jet-torus structure. Left: Crab Nebula. Center: Vela. Right: B1509 



SN ejecta. The SN ejecta are in free expansion, so this phase is referred as free 
expansion phase. During this phase the pulsar spin-down luminosity is high 
and almost constant, and the nebula can be easily observed in high energy 
emission. The Crab Nebula is presently supposed to be in this phase. The free- 
expansion phase terminates after about 1000-3000 yr, when the PWN reaches 
the reverse shock in the SNR shell. From this moment, the evolution of the 
PWN is deeply modified by interaction with the more massive and energetic 
SNR ejecta. In the absence of a central PWN, the reverse shock is supposed 
to reach the center of the SNR in about 5000-10000 yr. As the SNR ejecta 
compress back the PWN, its pressure increases. The compressed PWN will 
then push back the inner edge of the SNR. This phase is called "reverberation 
phase" : the PWN might undergo several compressions and rarefactions. Even 
if the energy injection from the pulsar is negligible at these later times, PWNe 
can still be observed, due to the re-energization during compressions. Multidi- 
mensional studies have shown t hat the SNR-PWN in terface is Rayleigh- Taylor 
unstable during compressions (IBlondin et al.l . 1200 ll ). This can cause efficient 
mixing of the pulsar wind material with the SNR, and ultimately with the 
ISM. The third phase corresponds to the Sedov phase of the SNR: the PWN 
expands adiabatically inside the heated SNR. Given the absence of energy 
injection, a PWN in this stage is probably only observable as a faint extended 
radio source. Recent studies have shown that this simple evolutionary picture 
can be further complicated if one takes i nto account the pulsar kick veloc 



ity with respect to the ce nter of the SNR (Ivan der Swaluw et al.l . 120031 . 12004 



Bucciantini et all l2005al ). If the kick is st rong enough, PWNe can be revealed 



also at later time as bow-shock nebulae (IBucciantini et al.l . l2005al ). when the 



pulsar emerges from the SNR, and interacts directly with the ISM. 



The interest of the scientific community for PWNe has received a consider- 
able impulse from recent optical and X-ray images from HST, CHANDRA 
and XMM-Newton. These new data show that the inner region of PWNe 
is characterized by a complex axisymmetric structure, g enerally referred as 



jet-torus structure (F ig. 1). First observed in the Crab dHester et al.l. I l995 



Weisskopf et~aTll2000h . it has been detected in many other PWNe ( jPavlov et al.l 
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2003 
2004 



Gaensler et all 12002: iLu et all 12002; IRomani & Nd . 120031: ISlane et al. 



Camilo et all . I2004J ; IRomani et all . 120051 ) . It consists of a main emission 



torus, in what is thought to be the equatorial plane of the pulsar rotation, often 
associated with multiple arcs or rings, a central knot, located in the vicinity of 
the pulsar, and one or two opposite jets along the polar axis, apparently orig- 
inating very close to the location of the pulsar itself. It was immediately clear 
that the simple ID radial model by KC84 was completely unable to explain 
the observed features. Even if qualitatively the existence of a main torus (the 
first feature of the jet-torus structure to be detected in the Crab Nebula) could 
be explained by assuming that the energy injectio n, or the particles accelera- 
tion ef ficiency were higher in the equatorial plane (IBogovalov fc Khangoulian . 
2002al ). it wa s not possible to repr oduce the observed luminosity by using the 
KC84 model. IShibata et all (120031 ) was the first to realize that, to explain the 
difference in brightness between the front and back sides of the torus in the 
Crab Nebula, the post-shock flow velocity should have been ~ 0.3 — 0.4c, much 
higher than typical values for subsonic expanding flows. As the quality of data 
increased the inconsistency of the old model was made more clear. The exis- 
tence of an inner ring, separated from the torus by an under-luminous region, 
was incompatible with the assumption of a smooth flow from the TS to the 
edge of the nebula. To justify the presence of the luminous X-ray knot inside 
the under-luminous supersonic wind region, dissipation of the magnetic field 
was thought to be at work, as a way to accelerate high energy particle in the 
supersonic wind. The presence of a jet was perhaps the most int eresting feature 
for the theoretical issues it raised (ILyubarsky fc Eichlerl . l200ll ) . The fact that 
the jet originated close to the pulsar, su ggested the poss i bility of collimation 



of th e pulsar wind. H owever, theoretical (IBegelman &: Lil. Il994t iBeskin et al 



2005 



1998) and numerical (IContopoulos et all Il999l; IBogovalovl . l2001al ; iGruzinov 



Komissarov . 2003 " Bucciantini et al. . 2006a ) studies of relativistic winds 



from pulsars, failed to achieve any form of collimated energ etic outflow. Fi- 



nally photon index maps of Crab Nebula (IMori et all 120041 ). show that the 



spectrum flattens moving from the inner ring toward the main torus, while the 
standard model predicts a steepening due to synchrotron losses. Interestingly 
the symmetry axis of the jet-torus appears to correspond to the major axis 
of the nebula, suggesting that the toroidal magnetic field is a key element in 
shaping the inner flow. 



Only recently a theoretical understanding of a possible mechanism to explain 
the jet-torus structure has been presented, and it was soon recognized that 
magnetization and energy distribution in the pulsar wind were key elements. 
While the old KC84 model assumed an isotropic energ y flux in the wi nd, it was 
known for a long time that the asympto tic solution by iMiche] fll973h. re cently 



2006 



confirmed with numerica l simulations (IBogovalovl . l2001al ; iKomissarov 
Bucciantini et all l2006al ) , predicted a higher equatorial flux. An axisymmetric 
energy flux en hanced at the equator, produces an oblate TS with a cusp in the 



polar region (IBogovalov fc Khangoulian! . l2002al ; IBogovalov fc Khangoulyanl . 
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2002bl ). Given that the TS is oblique to the wind direction at higher latitudes, 
the post shock flow in the nebula can have speeds ~ 0.3 — 0.5c. Moreover if 
hoop-stresses are at work in the mildly rela t ivisti c flow, jet collimation might 
occur in the post shock region (ILyubarskyl . |2002| ; iKhangoulian fc Bogovalovl . 
20031 ). The evident complexity of this scenario clearly prevented any sophisti- 



cated theoretical model, however, the recent de velopment of more efficient and 



robust numerical schemes fo r relativistic MHD flKomissarovl . ll999l ; lDel Zanna et al. 



20031 ; iGammie et all 120031 ). has open the possibility for a detailed numeri- 



cal description. The initial numerical results (jKomissarov fc Lyubarskyl . 120041 ) 
show that, based on the theoretical assumption of axisymmetric energy flux 
and post-shock hoop stresses, it is indeed possible to reproduce the observed 
features. Subsequent studies have tried to refine this picture, by investigating 
if and how the emission features could be used to understand the conditions 
of the pre-shock wind, trying to constrain the nebular emission properties, 
polarization and spectral variations. Despite the success in explaining the for- 
mation of the jet-torus structure, the recent developments have also posed 
other questions, regarding the fine details and time variability, which have not 
yet been satisfactory answered. 



2 Global structure of PWNe 



Let us first consider the global structure and morphology of PWNe, in the 
first phase of the PWN-SNR evolution, when the pulsar energy injection is 
still high and the PWN expands inside the free-expanding ejecta of the SNR. 
This corresponds to the present phase of the Crab Nebula, and in general is the 
phase in which X-ray emission from PWNe is expected to be stronger. In Fig. 2 
we present a simplified picture of a spherically symmetric P WN together with 



the re sult of a relativistic hydrodynamical simulation by iBucciantini et al 



( 120031 ). Around the pulsar there is a region occupied be the ultra-relativistic 
pulsar wind. The wind is terminated in a strong shock. Downstream of the 
shock, there is a region filled with a low density hot plasma. The observed non- 
thermal emission comes from this region, which constitutes the PWN itself. As 
this region expands inside the high density, cold, supersonic ejecta of the SNR, 
a thin shell of compressed material is formed. The system is confined inside 
the SNR shell where the reverse shock, the forward shock, and the contact 
discontinuity, separating matter of the SN from the shocked ISM, are visible. 
In the following I will consider only the case of a pulsar at rest in the center 
of the system. The expansion velocity of PWNe in this early phase is known 
to be few thousands kilometer per second, while typical pulsar velocities are 
in the range 50-300 km/s. The pulsar kick will not appreciably modify the 
evolution and properties of the PWN before the onset of the reverberation 
stage. However, after the onset of the reverberation phase, major deviations in 
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Fig. 2. Left picture: schematic representation of the global structure of the PWN 
in the first phase of its evolution inside a SNR. Right picture: density profile ob- 
tained from a time dependent numer ical simulation in the relativistic hydrodynam- 
ical regime ( Bucciantini et al. . 20031 ). From the center the various regions are: 1- 
the relativistic pulsar wind, 2- the hot magnetized bubble responsible for the non 
thermal emission, 3- the free expanding ejecta of the SNR, 4- the ISM. These regions 
are separated by discontinuities: a- the wind termination shock, b- the contact dis- 
continuity between the hot shocked pulsar material and the swept-up SNR ejecta, 
c-the front shock of the thin shell expanding into the ejecta, d- the reverse shock 
of the SNR, e- the contact discontinuity separating the ejecta material from the 
compressed ISM, f- the forward SNR shock. 



the global morphology and evolution of the PWN-SNR system can be expected 
(for a mor e complete study of the effect s of pulsar velocity on the SNR-PWN 
system see van der Swaluw et al. f 20041 )). 



Given that the edge of the PWN is bounded by a compressed thin shell 
it is possible to derive analytically the evolution of the PWN radius with 
time, if one assume the SNR ejecta to be in a self-similar configuration. It 
turns out, in this case, that the evolution of the PWN is self-similar too, de- 
pending only on the pulsar luminosity, the density and vel ocity of the ejecta 



( iReynolds fc Chevalier! . 1 19841 : Ivan der Swaluw et all l200ll ). Even if analytic 



models assume a uniform pressure in the PW N, an assumption valid o nly in 



the pure hydrodynamical regime, simulations (jBucciantini et all 120031 ) have 



shown, in the ID radial case, that the boundary pressure at the edge of the 
PWN does not depend on the magnetization, but only on the total energy 
injected by the pulsar. This suggests that global evolution of the size of PWN 
is not sensitive to the overall magnetization, thus strengthening the reliability 
of pressure driven thin shell models (Gelfand et al. these proceedings). The 
effect of spin-down and th e drop in pulsar luminosity have also been studied 
( jBucciantini et all l2004al ). Again, results show that the global evolution can 
be described analytically, using the thin shell approximation, and considering 
only the total energy injected by the pulsar. Simulations suggest however that 
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a strong spin-down might change the ratio between the radius of the TS and 
that of the nebula. Given that such ratio has been used in the past to infer 
the magnetization of the wind, results show that temporal evolution in the 
pulsar properties can significantly alter the inferred values. 



The results of both analytic theory and spherically symmetric numerical sim- 
ulations, show that the thin shell surrounding the PWN accelerates. Given 
that the density of the shell is much higher than the inertia of the relativistic 
plasma, the shell will be subject to Rayleigh- Taylor instability. The filamen- 
tary network of the Crab Nebu la, observable in lin e emission, is supposed to 



be the result of such inst abili t y ([Hester et al.l.ll996T). Various simulations have 
investigated its growth (IJunl . Il998l ; iBucciantini et all l2004j ). showing that if 



the field is highl y ordered at the bounda ry it can easily suppress the growth of 



parallel modes (IBucciantini et all I2004I ). Recent results (Stone, private com- 



munication) suggest that a more complex geometry of the magnetic field in 
the outer part of the nebula is required to reproduce the observed filamentary 
shape. 



Magneti zation plays a key role also in shaping the nebula, as s hown in the early 
work b y lBegefman fc Lil (Il992l ). and confirmed numerically by Ivan der Swaluw 
( 120031 ). Given that the magnetic field in the pulsar wind is mostly toroidal, if 
no instability arises to destroy the axisymmetric structure ( iBegelmanl . Il998l ). 
the magnetic field will pile up at the edge of the nebula. Due to magnetic stress 
the pressure driving the expansion of the thin shell is higher along the rotation 
axis of the pulsar than at the equator. Despite the change in boundary pres- 
sure, results show that the simple power law expansion rate, derived in the 
self-similar spherical model still holds. The elongation of the nebula increases 
with increasing magnetization, so that in principle can be used to constrain 
properties of the wind. Time dependent results show also that the elongation 
tends to saturate in a few hundreds years, a period smaller than the age of ob- 
serve d nebulae, and that the saturation time is sho rter for smaller magnetiza- 



tion (Ivan der Swaluwl. 12 003; De l Zanna et all 12004 ) . Interesti ngly, simulations 



done by iDel Zanna et al.l (120041 ) give results consistent with Ivan der Swaluw 
( 120031 ). despite the fact that in the former the pulsar wind is highly anisotropic, 
and the energy flux is concentrated in the equatorial plane. The reason is that 
even if the wind is anisotropic, and the post-shock flow has high velocities, 
in the bulk of the nebula the flow is subsonic and pressure equilibrium holds. 
This confirms that while wind anisotropy might be important in the model- 
ing of the inner X-ray emitting region, the global evolution and large scale 
properties of the nebula depend mostly on the integrated pulsar luminosity. 
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3 The inner structure of PWNe 



Let us now turn our attention to the inner region of PWNe, where the X- 
ray emission comes from. Simplified ID models, assume an isotropic energy 
flux, however, a series of recent numerical s imulations both in the force- 
free (Contopoulos et all Il999 : Gruzinov . 2005) an d in the MHD ( Bogovalov . 



2001al ; iKomissarovl . 120061 ; iBucciantini et all l2006al ) regimes have shown that, 



far away from the light cylinder, th e flow of an a ligned rotator closely resem- 
bles the exact monopole solution by iMichell (119731 ). This solution predicts that 
at the typical distance of the TS, the energy flux in the wind has a strong lati- 
tudinal dependence of the form L{6) = L Q (1 +a sin 2 (9)), where a is a measure 
of the pole-equator anisotropy (a » 1), while the magnetic field in the wind 
B{6) oc sin(#). In this case the magnetization in the wind is a function of 
the polar angle cr(9). Various numerical simulation s of the interaction of such 
wind with the SNR electa have been presen t ed (IKomissarov Lvubarskyl . 



2004 ; iDel Zanna et all 12004 ; iBogovalov et all l2005t iDel Zanna et all 120061 ) 



Interestingly these results show that, as long as the energy flux is the same, 
the post shock flow does not depend on the specific values of Lorentz fac- 
tor or density distribution. In this sense flow dynamics cannot be used to 
constrain the value of the wind Lorentz factor or composition. In Fig. 3 the 
shape of the TS is shown. Due to the higher equatorial energy flux the TS 
is closer to the pulsar along the polar axis than on the equator. In the low 
magnetization case (a eguator < 0.01) the shape of the shock can be reasonably 
reproduced by assuming pressure equilibrium between the wind and the neb- 
ula. However, due to magnetic stresses, the nebular pressure is higher on the 
axis, and the shape of the termination shock becomes more oblate at higher 
a. Due to the TS shape, the flow is slowed down to speed ~ c/3 close to the 
equator, but at higher latitudes, where the shock is oblique, the post shock 
flow is still super-fastmagnetosonic. Moreover, plasma at higher latitudes is 
deflected toward the equatorial plane. This flow is partly slowed in the so 
called "rim-shock", however its velocities are still supersonic. As shown in 
Fig. 3 the result of the anisotropic energy distribution is that almost all of the 
downstream plasma is confined in the equatorial region. The inner part moves 
with speed ~ c/3, while the outer channels has velocity ~ 0.5c, in agreement 
with the value inferred from the differ ence in the lum i nosity of front and back 



side of the torus in the Crab Nebula (jShibata et all 120031 ). Numerical simu- 
lations seem to give contrasting results regarding the stability o f the TS. In 



some cases the shape o f the TS evolves in a self similar manner (jVolpil . 12006 



Del Zanna et all 12004) . while in other cases it shows significant fluctuations 



( IKomissarov fc Lvubarskyl . 12004 ). Perhaps such behavior might be associated 
with boundary conditions on the expansion velocity and the average magne- 
tization in the nebula, but no detailed study has yet been presented on how 
the boundary velocity affects the overall inner flow pattern. 
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Fig. 3. Left panel: g lobal structure of a PWN, pressure distribution (from 
Del Zanna et al. ( 20041 )). Despite the higher equatorial energy flux, the pressure 



on large scales depends on the cylindrical radius as expected from pressure equilib 
rium. The higher polar pressure is responsible for the elongation of the nebula. Right 



panel : structure of the flow at the termination shock, velocity (from lDel Zanna et al 



The various labels indicates: A- the relativistic pulsar wind, B- the subsonic 



equatorial outflow, C- supersonic equatorial funnel, D- super-fastmagnetosonic flow, 
a- termination shock, b- rim shock, c- fastmagnetosonic surface. 

As the flow expands in this equatorial sheet away from the TS, the magneti- 
zation increases, until equipartition is reached. Given that the magnetization 
in the wind is higher in the equatorial region, the flow reaches equipartition 
close to the equator, before that at higher latitudes. The post-shock flow is 
only mildly relativistic, so hoop stresses can produce collimation. Given the 
extra degree of freedom with respect to the simple radial model by KC84 (the 
plasma can now move also in then 9 direction), once equipartition is reached, 
the magnetic pressure prevents further compression of the magnetic field, and 
the flow is diverted back toward the axis. This is the process that causes colli- 
mation and the formation of a jet along the axis itself. Fig. 4 shows numerical 
results for the magnetization and velocity profile inside a nebula. Obviously 
the values of the wind magnetization play a key role in the formation of the 
jet. If magnetization is too small, then equipartition is not reached inside 
the nebula, hoop stresses never become effective, the equatorial channel sur- 
vives to the edge of the PWN and no jet is formed. At higher magnetization 
equipartition is reached in the close vicinity of the TS, and most of the plasma 
is diverted and collimated in a jet. The size and flow velocity in the jet are 
a function of a. For values o ~ 0.03 the plasma speed in the jet is ~ 0.7c in 



agreement w i th ob s ervation of the jet in Crab Nebula (jWeisskopf et all 12000 



Hester et all 120021 ; iMelatos et all 120051 ). Simulations also show that there is 



a global meridional circulation inside the nebula associated with the jet for- 
mation , with typical speeds ~ 0.1c. Even if this speed is quite subsonic (as 
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Fig. 4. Result of a numerical simulation s howing the details of the formation of 
the polar jet (from iDel Zanna et al.l (|2004l )). Left panel: magnetization inside the 
nebula. Right panel: flow velocity. Notice that when equipartition is reached the 
flow is diverted and collimated toward the axis. 

previously mentioned in the context of global pressure distribution), it might 
have important consequences for the interaction with the SNR swept-up shell, 
and especially on the develo pment of l ocal s hear instabilities and on the pos- 
sible mixing with cold ions (iLyutikovl . 120031 ) . 



The picture presented above can be modified to take into account the fact that 
a puls ar, most probably, is an oblique rotat or. It has been sho wn both theoret- 
ically (jBogovalovl . l2001bl ) and numerically (jSpitkovskyl . 120061 ) , that even in the 
case of an oblique rotator, the energy distribution in the wind is essentially the 
same as the aligned case (in the case of split monopole it can be shown to be 
identical). However, the presence of a folded current sheet that might extend 
to higher latitudes, has important consequences on the magnetization. If the 
alternating magnetic field in t his striped wind region is dissipated, and this 
can h appen either in the wind (jLyubarskv fc Kirkl. l200ll; iKirk fc Skiaeraasenl . 



2003 



2005h . 



then the 



20031 ) or at the termination shock itself (jLyubarskyl . 
magnetization will drop to zero on the equator and will reach a maximum at 
intermediate latitudes, depending on the obliquity of the pulsar. Given that, in 
MHD codes, the TS is a thin discontinuity, there is no difference between dissi- 
pation in the wind or at the shock, and one can chose the former for numerical 
simplicity. We need however to remember that, in simulations of aligned ro- 
tators, the basic assumtpio is that the equatorial current sheet survives in 
the nebula. In reality in the shocked plasma the current sheet can undergo 
instability and dissipation, the details of which have not being investigated. 

The presence of an unmagnetized region close to the equator, adds complex- 
ity to the picture presented above. In particular, now the equatorial flow does 
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Fig. 5. Effects of a striped wind region: result of a numerical simul a tion s howing 
the detail of the formation of the polar jet (from Del Zanna et al. ( 20041 )). Left 
panel: magnetization inside the nebula. Right panel: flow velocity. An equatorial 
fast outflow remains, while only material at intermediate latitudes is collimated 
toward the axis. 

not reach equipartition inside the nebula, and an equatorial channel survives 
to much larger distances from the termination shock. However, at higher lati- 
tudes, where dissipation is marginal, and the magnetization is higher, equipar- 
tition can be reached and the flow is then diverted to the axis and collimated 
into a jet. Fig. 5 show the flow structure. As will be discussed in more detail 
in the following, emission maps based on the results of numerical simulations 
suggest that a striped wind model is more promising in explaining the ob- 
serves inner-ring outer-torus structure of many PWNe. The extra complexity 
added by the presence of a striped wind region, can b e used as a probe f or the 



properties of the wind. In particular recent results (IDel Zanna et all 120061 ) 



have shown that, given the same average magnetization, a bigger striped wind 
region, leads to weaker jets, less intense hoop-stresses, and a bigger size of the 
termination shock, all of which are observable. 



Even if the fluid structure recovered in relativistic MHD simulations clearly 
shows the formation of a jet, to properly assess their validity one must consider 
the emission and compare results based on simulations with observations. De- 
spite the fact that, variations in the particles' energy distribution and accelera- 
tion along the TS are poorely known, to a first order approximation one might 
assume a power-law injection distribution, and an uniform efficiency. Under 
this assumption, to obtain X-ray maps, one must then com pute the adiabatic 
and synchrotron losses. In iKomissarov &: Lyubarskyl (120041 ) this was done by 
assuming an exponentially decaying factor, with spatial constant cooling time 
depending on the frequency of observation and an average magn etic field. A 
more sophisticated strategy was adopted by IShibata et al.l (120031 ) where sev- 
eral particles energies are evolved independently along streamlines. While the 
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Fig. 6. Emission f eatures. Left panel: optical HST image of the inner region in 



the Crab Nebula (|Hester et all Il995l ). Central panel: optical synchrotron map 
based on a numerical simulation of the flow inside a PWN, in the case of a 
strip ed wind. Schematic interp r etatio n of the origin of some of the observed fea- 
ture ( Komissarov &: Lyubarsky . 20041 ) . The various labels indicate: 1- the knot, 2- 
the wisps, 3- the main torus, 4- the anvil, 5- the back side of the wisps. Notice that 
this features are recovered in the simulated map. 

latter is probably the more accurate, it can be applied only in steady state 
cases, where the streamlines ca n be traced without am biguities. An interme- 
diate approach was proposed by lDel Zanna et al.l (120061 ). where the maximum 
energy e max of the particles' distribution was evolved with the flow, accounting 
both for synchrotron and adiabatic losses. In this case the particle distribution 
remains a power low with a break at e max . 

Preliminary results show that it is possible to reproduce both the inner ring (or 
wisps region) and the outer torus, only by assuming the presence of a striped 
wind region. Indeed in the case of no striped wind, maps usually produce a 
single luminous arch. In Fig. 6, we show an optical image of the Crab Nebula, 
a map based on a simulation with striped wind (in this case, for comparison 
in th e visible band, synchrotron coolin g is neglected), and a schematic picture 
from iKomissarov fc Lyubarskyl (120041 1 , to illustrate the origins of the various 
features. The knot and the inner ring are both explained as due to the high 
velocity flow in the immediate post shock region, at intermediate latitudes. 
However, by looking at Fig. 5, we see that, in the striped wind fast 
equatorial channel survives, extending far from the TS. It is the emission com- 
ing from this channel that is respo nsible for the main torus observed in the 



Crab Nebula. Recent investigations (jVolpil . 120061 ; iDel Zanna et all 120061 ) have 



shown that the size and shape of the striped wind region, have important ob- 
servational consequences. In general a larger unmagnetized sector will produce 
a more structured inner emission; the region corresponding to the inner ring 
might split in several minor arches; the extent of the X-ray nebula is bigger; 
the jet is thinner and less enhanced. If the magnetization is not high enough, 
energetic particles can survive all along the jet and accumulate in a mushroom 
cap at the edge of the nebula. Simulations also seem to disfav or configuration 
with a higher magnetization close to the pole (lAronsl . Il998l ). In fact in this 
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runA t = laCOyr o=D.B : Photon Index 




Fig. 7. Spectral p roperties of PWNe . Left panel: map of the X-ray photon index in 
the Crab nebula (jMori et all l2004h . Right panel: photon index computed from a 
numerical simulation. Dissipation was added in the polar region to reproduce the 
correct value for the jet. 

case, hoop stresses are less effective and no observable bright extended jet 
is formed. While modeling of the equatorial region does not need any accu- 
rate treatment of the sy nchrotron losses (the feature are observed also in radio 
( iBietenholz et al.L I2004T )). these are fundamental in understanding the jet, and 
simulated map clearly show that the jet is not visible over the background, in 
optical or radio. 



Interestingly, a detailed study of emission with spectral maps, polarization 
properties, and integrated spectra, can also be used to better constrain the 
existing models. If the particle energy distribution is somehow followed along 
a streamline, it is possible to reconstruct spectral maps. In Fig. 7 an X-ray 
map of the spectral index based on simulati ons is presented together with a 
map of the Crab Nebula (jMori et all 12004 ). The main features are indeed 
recovered. We also see that simulations produce maps where the spectrum 
appears to flatten moving toward the main torus, without the need to as- 
sume any re-acceleration. This is just due to the Doppler boosting effect. At 
low speeds, when observing at a given frequency, we are implicitly sampling 
particles all at the same energy (let us neglect variation of magnetic field). 
When velocities are relativistic, the particle energy responsible for emission at 
a given frequency, depends also on the Doppler boosting. The spectrum ap- 
pears to flatten in the torus, because, given the higher speed, we are looking 
a lower energy particles. However, map based on simulations fail to recover 
the correct spectrum in the jet: the simulated spectrum is much steeper than 
what is observed, due to excessive synchrotron losses. To reproduce the obser- 
vations, one must assume some form of dissipation and re-energization along 
the axis. This could be associated with local instabilities of the toroidal mag- 
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Fig. 8. Effect of th e striped wind on the X-ray emission features (from 
Del Zanna et al.1 (|200d )). Both maps are based on simulations with the same aver- 
age wind magnetization, but the one on the left has a much smaller striped wind 
region. Notice the different extent of the emission, the presence of more arches in 
the right image, and the luminosity of the jet. 



netic field ( iBegelmanl . Il998l ) which are not captured by a xisymmetric s i mula- 



tions, but for which there is much observational evidence (jPavlov et al.1 . 12003 



Mori et all 12004! : iMelatos et all . 120051 : IPeLanev et~aD . |2006|). By considering 



the integrated spectra we can obtain some constraints on the magnetization. 
Preliminary results show that o > 0.03 is required and a configuration with a 
smaller striped wind region are favored. Moreover, integrated spectra show a 
flattening above the CHANDRA band, which has been interpreted as due to 
the high velocity in the post-shock region. In this case in fact, the size of the 
nebula at higher energies, does not seem to drop substantially, as it does from 
optical to X-ray. 



Polarization might also be used to derive useful information about the mag- 
netic field structure. Emission maps depend strongly on the flow velocity inside 
the nebula, but they are not very sensitive to the presence of a small scale dis- 
ordered component of the magnetic field. It is thus not possible at the moment 
to use them to properly constrain the geometry of the field. On the contrary, 
polarization might prove very useful. Given that only optical polarization is 
available, any study should be limited to the brightest features, like the wisp, 
inner ring and tor us. The effect of fl o w velo city o n the polarizatio n angle has 
been discussed by iBucciantini et al.l ( j2005bl ) and iDel Zanna et al. 16x3). By 
comparison with old optical polarizatio n maps, one can imme diately see that 
some of the properties are recovered. In ISchmidt et al.l (119791 ) it is also inter- 
esting to note the presence of a region which can be identified with the jet 
of Crab. More recent data (Graham, private communication) have also shown 
that the knot has a high degree of polarization, and that the polarization angle 
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is consistent with this feature originating in the post shock flow. 



4 Time variability 



Most of the work done in the past has focused on the explanation of the 
main observed features that appear to be quite persistent on long timescale. 
However, there is increasing interest in the understanding of time variabil- 
ity. Variability of the wisps in the Crab Nebula has been known since their 
discovery (see refer e nces in lAronsI (119981)'). and high resolution observations 



( jHester et all 12002 ; iBietenholz et all 120041 ) have allowed us to stu dy its be- 



havio ur in detail. The jet in Vela also ap pears to be strongly variable (IPaylov et al. 
20031 ). as do es the jet of both the C rab (Mori et all 120041 ; iMelatos et all l20olh ~ 



and B1509 (lDeLaney et al.l . 120061 ). B1509 shows also variability in the inner 
ring. There is also evidence for the presence of variable knots (non axisym- 
metric features) both in Crab and B1509. Variability in the jet, which usually 
has timescale of years, has been associated with kink or sa usage modes in 



the s trong toroidal field, or even with the fire-hose instability (ITrussoni et al. 



19881 ) in the case of Vela. Recently in B1509, yariable knots have be en observed 
inside the TS, on the side of the counter-jet (lDeLaney et al.l . 120061 ) . These are 
non axisymmetric features, however it has been speculated that they might 
be the signature of instability in the TS cusps at the base of the counter-jet. 



More interesting is the case of the wisps, which show variability on a timescale 
of weeks. The presence of high velocity flow channel s inside the nebu la sug- 
gests the possibility of Kelvin-Helmholtz instability (iBegelmanl . Il999l ). how- 
ever it seems that it cannot account for the lumino sity variations required in 
the wis ps (the fluctuations seem not to exceed 10%. iBucciantini fc Del Zanna 
(l2006bh ). Towever, the flow inside the nebula is strongly dynamic, and varia- 
tions in the meri dional circulation patter n, can provide a possible explanation. 
In simulations by lBogovalov et al.l (120051 ) th e synchrotron luminosity varies on 
short time scales (few months). Results by iKomissarov & Lyubarskyj (120041 ) 
also suggest that the flow might have a feedback action on the TS, causing it 
to change shape, thus inducing a change in the wisps, which originate in the 
post-shock flow. 



A diffe rent model for the variability of the wisps has been proposed by lSpitkovsky fe Arons 
( 120041 ) . The model is based on the assumption that high energy ions are 
present in the equatorial region of the wind. There are several reasons for 
this claim, most important the fact that, in kinetic simulations of acceleration 
in a str ong shock, a pure pair plasma does not produce a power-law distri- 
bution (lAmato fc Aronsl . 120061 ) . Given that ions have a much larger larmor 
radius, of order of the size of the wisp region, they introduce a substantial 
deviation from a pure fluid picture. In particular, electrons are compressed by 
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the gyration of ions and the emission is enhanced. The model is able to repro- 
duce the observed time scale variability and the average distance between the 
torus and inner wisps. However, these results are based on a ID spherically 
symmetric model, do not take into account the energy flux anisotropy, and 
the shape of the TS, and assume coherent motion of the ions in an equatorial 
sector. The ion cyclotron and MHD model are by no mean in compatible. If 



the presence of high ene rgy ions is confirmed by observations (lAmato et al. 



20031 ; iHorns et all 120061 ) . an hybrid approach that integrates both features 



would be necessary fo r better model s . A th ird model for the wisps dynamics 
has been presented by iFoy fe Hester! (120061 ) , based on the idea of synchrotron 
cooling instability. However, in order to work, this model assumes most of the 
energy in the nebula to be carried by the X-ray emitting particles, in contrast 
with the standard assumption that lower energy particles dominated the bulk 
flow. 



5 Conclusion 



In the last few years, thanks to a combination of high resolution observations, 
and numerical simulation, our understanding of the evolution and internal 
dynamics of PWNe has greatly improved. We have been able to reproduce 
the observed jet-torus structure and to clarify how the jet is formed in the 
post shock flow, as a function of the wind magnetization. Simulated maps can 
reproduce many of the observed features, and the correct behavior of the spec- 
tral index in X-ray. Results suggest that, the best agreement is achieved in the 
case of a striped wind, even if MHD simulations are not able to distinguish 
between dissipation of the current sheet in the wind or at the TS. It appears 
that the size of the striped wind zone affects both the complexity of the ob- 
served structure and the intensity of the jet. This has open the possibility of 
using X-ray imaging to constrain the pulsar wind properties; already the rings 
and tori obser ved in many PWNe have been used to determine the spin axis 



of the pulsar (IRomani et all 120051 ) . Interestingly in the Crab the inner ring 
does not appear boosted, while the wisps (which are interpreted as its optical 
counterpart) are. 

There are still however unsolved questions, and possible future developments 
for research in this field. The interaction of the relativistic hot plasma with the 
confining SNR, has not been considered in detail, in the framework of the new 
interior models. It is not clear how the Rayleigh- Taylor instability is affected by 
the global circulation, and what kind of feedback is exerted by the filamentary 
network on the inner flow. All simulations are axisymmetric, thus none is able 
to address the problem of the stability of the toroidal field. Observations on the 
other hand show clear non axisymmetric variability in the jet. Such instability 
might be important in modeling the jet emission, especially given that present 
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results fail to reproduce the correct synchrotron spectrum. It is not clear if 
small scale disordered field is present in the inner region (may be a residual of 
the dissipation in the TS of the striped wind). A combination of simulations 
and polarimetry might help solving this question. Finally it should be borne 
in mind the possible presence of ions and the modifications they will produce. 
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